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ABSTRACT 


Pe study of the ae and second-order, dynamic, ideal 
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are obtained with a FORTRAN IV computer program based on the 
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First-order added mass and damping in heave for a semi- 
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are presented. Second-order solutions for the semi-circular 
Meeeorin in eave also including finite depth effects are 
obtained. The solutions are verified by comparing with 


existing theory and experiment. 





TABLE OF CONTENTS 


as LM ERODIC{TON (6 25 16 
A. HISTORICAL BACKGROUND ---------------------- 16 
1. First-Order Solutions ------------------ uly 
a. General Potential Theory ----------- Li 
b. Finite Element Solutions ----------- i 
2. Second-Order Solutions ----------------- 20 
a. General Potential Theory ----------- 20 
b. Finite Element Solutions ----------- Ze 
B, OBJECTIVES OF PRESENT RESEARCH ------------- 25 
a. FORMULATION OF THE MATHEMATICAL MODEL ---------- oy 
A. THE NONLINEAR BOUNDARY VALUE PROBLEM ------- oy 
B. PERTURBATION ANALYSIS ---------------------- 29 
1. The Linearized Free Surface Condition -- 30 
2. The Ship Interface Condition ----------- Sl 
Bee rerturbarion Formulas for Pressure, 
Force, and Wave Amplitude -------------- 35 
C. THE BOUNDARY VALUE For o‘1) ~--------------- 39 
D. THE BOUNDARY VALUE PROBLEM FoR 92) -------- 39 
BE. THE SECOND-ORDER BOUNDARY VALUE 
PROBLEM WITH FINITE DEPTH ------------------ 4 
1. The Second-Order Boundary Value 
Problem for 762) —------—---------------- 42 
F, FORMULAS FOR PRESSURE, FORCE, AND 
ee en ae 43 
III. FINITE ELEMENT DISCRETIZATION ------------------ 45 
A, THE DISCRETIZED FIRST-ORDER PROBLEM -------- 47 





B. THE DISCRETIZED SECOND-ORDER PROBLEM --------- 48 
1. Calculation of Nodal Values of the 
Field Derivatives of gS)? ee eee eee 50 
TV. THE COMPUTER PROGRAM AND NUMERICAL SOLUTION ------ 52 
A. MAIN FEATURES OF THE COMPUTER PROGRAM -------- 52 
Pewee cORUTHONS —---------------___.-..--- 5 
1. First-Order Numerical Solutions ---------- 54 
a. Mesh Requirements -------------------- 54 
Pe shmiaesebceparation =——-—<=--—---=.=—-= 56 
c. Convergence of First-Order 
Solutions ---------------------------- Sif 
2. Second-Order Numerical Solutions --------- Sy i 
a. Mesh Requirements -------------------- 51 
b. The Second-Order Boundary Condition 
on S23 =--~-~~--~~-~--- -----_--_---- ----- 58 
ec. The Ship-Fluid Interface Second- 
Sudemeleciiaaie condi a.on ———————-——_ —_—— 60 
dad. Convergence of Second-Order 
Salute one 2 = 61 
Vv. SURI, Rei = 62 
A, DEFINITIONS ---------- ~-—-----~-------~-------- 62 
1. Added Mass and Damping Coefficients ------ 62 
2. Dimensionless Force Coefficients --------- 64 
3. Wave Amplitude at Infinity --------------- 65 
2, UI, FORMS 6.55. 65 
C. FIRST-ORDER NUMERICAL SOLUTIONS -------------- 67 
1. Semi-Circular Hull ----------------------- oT 
®, Bulb Bui 22 73 
en ome mer awoemeo und? es.)————- - --_- —_-______ —— is 





PSEOCND-ORDER NUMHRICAL RESULTS -------------~- 17 


1. Semi-Circular Hull — Infinite Depth ----- 77 

2. oemi-Circular Hull — Finite Depth ------- fs, 
vi. CONCLUSIONS ------------------------------------- 88 
Pier Pemoonplex Algebra —----------.-.---.------_-- 90 
APPENDIX B. Isoparametric Elements ------------------- 91 
APPENDIX C. Computer Program Economization ----------- 95 
BIBLIOGRAPHY ------------------------------------------ 98 
INITIAL DISTRIBUTION LIST ----------------------------- ILO 





mast 1. 


meee Il. 


LIST OF TABLES 


DIMENSIONLESS FORCE COEFFICIENTS ---------- 


CONSTANTS FOR HULL FORMS 





Figure 


Oe 


Bel. 


Jha 


13. 
14, 


i 


BEG 


17. 


iit eos Hh LOURES 


Title 
Hull at Free Surface ------------------------- 
Reduced Region ------------------------------- 
Hull Coordinate System ----------------------- 
Finite Element Mesh -------------------------- 
Hull Forms ----------------------------------- 


Added Mass and Damping Coefficients in 
Heave, Semi-Circular Hull, Infinite Depth ---- 


Variation of Added Mass Coefficient with 
Depth, Semi-Circular Hull in Heave --~--------- 


Variation of Damping Coefficient with 
Depth, Semi-Circular Hull in Heave ----------- 


Added Mass and Damping Coefficients in Sway, 
Semi-Circular Hull, Infinite Depth ----------- 


Exciting Force Amplitude in Heave, Bulb 
Hull, Infinite Depth ------------------------- 


Wave Amplitude for Heave, Bulb Hull, 
Infinite Depth -----------------~------------- 


Variation of Added Mass and Damping 
Poem wercits win Depyn, Bulb Hull in Heave -- 


Reduced Region for Wave Generator Problem ---- 


Wave Amplitude Generated by a Plunger Type 
NAS GSMNZEII CR? 65 eerer re aaa 


Wave Amplitude Generated by a Flapper Type 
asge GSQIEEtOle 2 ee 


Second-Order Exciting Force Amplitude in 
Neave,octit—Ccirecnular Hult, Infinite Depth ---- 


Second-Order Exciting Force Phase Angle 
PMimlcave ro ocmm—-Cimcular Hull, Infinite Depth - 


80 


81 





faeure 
FO . 


io. 


20. 


2ge 


Ee. 


ie ie 


ol. 


ace Le 


Variation of Second-Order Exciting Force 
Amplitudes with Depth, Semi-Circular 


Mey dys, Hieeawea 666653533332 =e 


Variation of Second-Order Exciting Force 
Phase Angles with Depth, Semi-Circular 


UU she biaeuwe (3S 


Ratio of Second-Order Exciting Force 
pimp breude Go First-Order Exciting Force 


Amplitude, Semi-Circular Hull in Heave ------- 
Total Exciting Force vs. Time ---------------- 


EeW@enGeerrdere wave Norm Effects —-—j—-—-———-———-———-— 


PieonenusMapping from €,n Plane to x,y Plane 


Nanaia momotecoetticitent Matrix —~—-—-———-=———~——~— 


IO 





symbol 


©2 2 i 


Ct ) 
pL) p(2) 


(1) =(2) =(2) 
eer. 4k. 


Ba) (1) 
x Vy 


p(2) p(2) 


2) (2) 
He ote 


Pas Ol vor NBOLS 


De faseiae Lon 
One-half of ship's cross-sectional area 
Amplitude of ship's motion (sway or heave) 
Ship's half-beam measured at y = 0 


Added mass and damping coefficients, 
respectively 


Celerity of first-order gravity wave and 
Stokes second-order wave 


Celerity of second-order gravity wave 


from velocity aGhemetar oo! 
Sap Searart 


Real coefficients in Lewis hull form 
equation 


Total nonlinear ship's exciting force 


Dainvmonieamsceeme—Ooracr components of F, 
respectively 


Nondimensional first and second-order 
force coefficients 


Resultant first-order hydrodynamic force 
acting on the ship (sway or heave, 
respectively). 


First and second-order force amplitudes, 
respectively. 


Amplitude of second-order static and 
dynamic components, respectively 


Acceleration of gravity 
Frequency dependent coefficients 
Fluid region depth 


Unit vectors in cartesian coordinate 


11 





elec, 0) ,P(x,0) 


aly) 
qty) 


ap ly) 


Py oho 8 


mie, 0) 


Demand t Lon 
The imaginary unit v-l 


First and second-order wave numbers, 
respectively (k, = o/c,) 


One-half of the total arc-length of a 
Syimeur.c Shap hull form 


TOGalenumoer Of nodes in fluid region 
finite element mesh 


Added mass 
Unit outward normal vector from region R 


Hydrodynamic overpressure (referred to 
atmospheric) 


zero, first, and second-order components 
Ol p 


Asymptotic functional forms 


Function describing mode shape of wave 
generator 


Function describing mode shape of plunger 
type wave generator 


Higeutonmoacsicr bing mode Shape of flapper 
type wave generator 


Moduli and arguments of complex numbers 


Function representing second-order Stokes 
free surface boundary condition 


The domain of the boundary value problem 
Houidarmy som the Region RK 
Denotes the ith region boundary 


Arc-length along ship's hull measured from 
Ship's centerline 


Time 


tyecomooncmus Ol tluid velocity vector, 
respectively 


WZ 





Symbol 


W 
min 


X5VJ 


me) ,y(s) 


XH yh 


€ 
mx, t ) 


q (1) xc) nb) (x,t) 
A? 


Ail) 72) 


En 


Pein c 1 Or 
ohip's velocity vector 
Reduced region width 
Reduced region width parameter 


Two-dimensional inertial cartesian 
coordinates 


Pocameurrc mull coordinates 

Functions of time such that exp,eyz, 
represent ship's harmonic motion in sway 
Ameo er respecrlvely 

Complex form of Yas 

Complex numbers 

Reduced region width parameter 

First and second-order exciting force 
phase angles meaSured with respect to 


ship's "top dead center" position 


Dimensionless frequency parameter 
(6 = o°b/g) 


Complex mapping parameter in Lewis hull 
form eduarpion 


Perturbation parameter (€ = a/b) 


Wave amplitude measured from y = 0 


First and second-order components of n(x,t) 


Dimensionless first-order wave amplitude 
Sno. alicia peaben', 


First and second-order dimensionless wave 
amplitude coefficients 


(HF) = FO) (x,¢)) 
Maximum amplitude of wave generator 
Cartesian coordinates in&,n plane 


MaeeGl Gens ..t y 





Symbol 
0 


mex .y st ) 


o0x,y,t) 


p61) 9(2) 


ge (1) 952) 


v 
wo 
v 
2 &) 


Ce) 


UA 2G UG 


CZ 


Definition 
Camerelagrerrequency Of Ship's motion 
Real=<valued nonlinear velocity potential 


Real-valued nonlinear velocity potential 
of a traveling gravity wave 


First and second-order components of 
Ceereciaoe ce ui Ve Ly 


First and second-order components of 
®*, respectively 


First and second-order complex, time- 
tidependemh. velocity potentials, 
respectively 


Modified complex, second-order velocity 
potential 


A known test velocity potential 


Second-order, complex Stokes wave 
velocity potential (time-independent) 


Meeecoluni Vecuor OL nodal values of 


g(t) c2)) 


or 
m xX m real, symmetric coefficient matrices 
Jacobian matrix (2 x 2) 


Total coefficient matrix (m x m) of the 
discretized first-order system 


1 x m row vector of cubic isoparametric 
Shape functions 


MeVbeomrPOn vector Of cubic 1sSOoparametric 
shape functions at the element level 


-acelumnn Vectors 
Mie vecvoremeradient operator 


Laplacian operator 


14 . 





ACKNOWLEDGEMENTS 


The author wishes to express his deep appreciation to 
Deeenobert &. Newton, Professor of Mechanical Engineering, 
for his advice and perceptive guidance throughout the 
development of this research. His patience and insight were 
key factors in the successful conclusion of this work. 

The author is indebted to Professor C. J. Garrison and 
Professor T. Sarpkaya for their helpful discussions during 
m@meecourse Of this study. 

The W. R. Church Computer Center of the Naval Postgraduate 
school provided the excellent computer service and technical 
eevee nor the numerical solutions. 

amelivemeune author wishes to thank his wife, Wendy, for 
maa oaulemoec understanding, and encouragement during the time 


of this research. 


i) 





te LNT RODUCTION 


Naval architects and marine engineers have long been 
Beeapey Jncerested in the study of the dynamic fluid response 
to the forced oscillations of a two-dimensional floating body. 
The applications of the studies include studying the rigid- 
body response of the ship due to an incoming wave train, and 
determining the natural frequency and vibration mode shapes 
Seemsnsp's Structure in connection with propulsion systems 
Besign [34]. 

Marine engineers, in modern times, are concerned with, 
among other things, the natural frequencies and stability 
characteristics of floating platforms near a coastline. 
fweetee they are concerned with the related problem of finding 
the dynamic loads on submerged, moored, or fixed bodies [8]. 

The modern techniques for solving these types of problems 
involve seeking the solution of a steady-state, periodic, 
potential flow problem with a free surface ina fluid of 
eitner infinite or finite depth. This problem is made 


tractable by linearizing the free surface boundary condition. 


A. HISTORICAL BACKGROUND 

the historical development of potential flow problems 
involving a free surface with wave motion, and the related 
problem of the dynamic response of fluids due to the motion 
of rigid bodies on and under a free surface is lengthy; but 


of interest. Several authors have reviewed this history in 
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pememdevarl (3,17,2e4,31)}. The history is repeated in the 
mollowing in the interest of providing a foundation upon 
feet this author endeavors to add additional information in 
oleae leeIipy FO Drovide a more complete history. 
oreo order Solutions 
a. General Potential Theory 

The foundation of the modern day theory of fluid 
ferens Desins, of course, with the classical works of Newton, 
Lagrange, Euler, Stokes, and others. 

However, in 1879, Lamb [15] provided one of the 
merseecomplete Studies of the theory of wave motions; which 
he updated five times until 1932. It appears that the next 
major contribution was made by Havelock [10] in 1929; wherein, 
he studied the solution of the wave forms created by a wave 
eenerator. 

ie medern beginnings of the study of the inter- 
peclon ot a fluid and a harmonically oscillating rigid body 
at the free surface start with the work of Lewis [19] in 1929. 
toweoeovudied the determination of added mass with a boundary 
SendiviOn requiring zero pressure at the mean position of the 
free surface, thus eliminating the effects of frequency of 
Dody oscillation and of damping. Lewis considered the high 
frequency case of the more general probiem. 

The first work which included the effects of 
frequency and damping was done by Ursell [38] in 1949. This 
work solved the problem of a right-circular, semi-~immersed 


cylinder heaving on the free surface of a fluid infinite in 
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extent. Ursell presented added mass and damping coefficients 
as a function of frequency. His method of solution involved 
mocing an “infinite” number of singularities of all orders 
Suen Ttiversection Of the free surface at rest and the 
Syeimaer centerline. The strengths of the sources were 
determined from the kinematic boundary condition applied on 
the immersed surface of the cylinder. Grim [9] arrived at a 
Semeetar solution in 1953. 

A period of about ten years passed before any 
additional work was reported on the subject. There is a 
rather obvious reason for this time lag. Ursell's solution 
maolved laborious numerical calculations which had to be 
done by hand. The advent of the digital computer made it 
Feasible to pursue further solutions. 

In 1960, Tasai [36] extended the work of Ursell 
to include roll and sway motions of Lewis forms. At the same 
time, Porter [30] compared linearized pressure calculations 
With experimental results, and extended the range of solutions 
Mmemnclude other hull forms. . 

It appears that Yu and Ursell [42] provided the 
first theoretical study of the effect of finite depth on 
two-dimensional solutions in 1961. This was followed by 
additional work by Kim [14] in 1969, who also included finite 
Seoen. in the first work, added mass coefficients and the 
properties of the resultant linear wave were reported. The 
peecond work, by Kim, gave added mass and damping coefficients 


aS a function of frequency and depth. Work by Paulling and 
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Richardson [28] and Paulling and Porter [29] provided some 
experimental verification of the theory in 1962. Vugts [40], 
at the Netherlands Research Center in 1968, reported additional 
feuemicive, Experimental research on seven different ship hull 
ferns . 

b. Finite Element Solutions 

Zienkiewicz [46], in 1965, provided the first 
Successful application of the finite element method to Poisson's 
equation. Although the finite element method was originally 
developed to solve problems in the theory of elastic mediums, 
Zienkiewicz's work "opened the door" to possible solutions to 
Bemicary Value problems of many different types, including the 
Peiution of fluid flow problems; this follows since Laplace's 
equation is contained in Poisson's equation. 

In the same year (1965), Zienkiewicz, Irons, and 
Nath [44] first used the finite elenwnt method to find added 
mass. Just as in the case of Lewis and others before, surface 
Waves were excluded from this solution, as well as in subse- 
quent works by Réren [32], Holand [12], Matsuura and Kawakami 
[22], and Matsumoto [21]. 

Zienkiewicz and Newton [43] presented the first 
general theoretical treatment of the fluid-structure inter- 
action problem using the finite element theory in 1969. Their 
work provided a means of including surface waves and deter- 
mination of energy loss due to wave propagation. The latter 
was accomplished by providing a radiation boundary condition. 


This radiation boundary condition was essential to the 
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Successful modeling of infinite regions under steady-state 
mem@eapsons. iin addition, compressibility effects could be 
included in the analysis. 

Sicmoudue¢i,94m 19/0, appears to be the first to 
Peemeecone finite element analysis developed by dienkiewicz 
and Newton to find added mass and damping coefficients as a 
Miiiemmbenl Of frequency and hull form. The fluid regions he 
Studied were chosen to simulate infinite depth. Some of his 
results are published in Ref. [23]. Bai [3] also included 
surface waves in his 1972 analysis of a wide range of problems. 
Mmemwork demonstrated the flexibility of the method by studying, 
am addition to the problems of Chenault, the additions of sway 
and roll motions, finite depth, irregular fluid bottom geometry, 
axisymmetric geometry, and fluid stratification. Bai also 
Meereoed criteria for properly placing the radiation boundary 
and developed a new form of the radiation condition for the 
axi-symmetric case. 

igemaumnOmemas Cxbended the studies of Chenault 
and Bai in the solution of the first-order problems with an 
improved version of Chenault's computer program. The results 
Peempresentved later in this work along with results separately 
meperted by Newton, Chenault, and Smith [23]. 

fee -econd Order Solutions 
a. General Potential Theory 

The origins of the second-order theory of gravity 

waves begin again with the work of Stokes and his classic 


Solution of nonlinear traveling waves. Others have extended 
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m@emwork Of stokes, but it appears that a modern version 

of Stokes' theory was first compiled by Stoker [35] in 1957. 
The generally accepted terminology associated with Stoker's 
work is "perturbation theory" which of course has applications 
in many fields. 

Fontannet [5], in 1961, appears to be the first 
GO attempt to use the basic nonlinear perturbation theory as 
described by Stoker to solve the problem of finding the 
characteristics of waves generated by a plane wave-maker using 
Lagrangian coordinates. Ogilvie [25], in 1963, found the 
first and second-order forces on a fixed or oscillating 
circular cylinder submerged in a fluid with a free surface. 
mewever, his second-order forces were only those due to the 
first-order solution, and only the time-independent part was 
reported. Tuck [37] found the second-order forces on a 
submerged cylinder in a uniform stream in 1965. Salvesen [33] 
extended Tuck's work to include wing-shaped bodies in his 
Ph.D. dissertation in 1966. 

Lee [17] and Parissis [27] appear to be the first 
tO apply second-order perturbation theory to attempt to solve 
Lhe fluid-structure interaction potential flow problem with 
a free surface in 1967. Lee confined his studies to the 
heaving of floating circular and U-shaped sections, while 
Parissis studied only the circular shape in heave; both 
studies were in fluids of infinite extent. These investiga- 
tions, as well as those previously mentioned, required the 


development of a fluid-structure second-order boundary 
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condition (discussed in more detail later) in addition to 
the free-surface equations compiled by Stoker. 

The work of Lee and Parissis was followed by that 
of Potash [31] in 1970. Potash extended the second-order 
studies to include the additional degrees of freedom of roll 
and sway and the attendant coupling between modes of motion. 
Recently, Garrison [7] completed a second-order solution of 
the problem of determining the dynamic loads experienced by 
a "fixed" three-dimensional body due to an incoming wave 
train. His theory includes submerged and surface-piercing 
meaLles . 

All of the second-order solutions previously 
mentioned used some form of singularity distribution in the 
meow tield. The functional forms of the singular functions 
used were uSually taken from the work of Wehausen and Laitone 
C41]. 

Db. Finite Blement Solutions 

The author has been able to find only one 
meme On FO the Second-order, of the class of problems being 
mmeeeuocea, Dy the finite element method. Allouard and Coudert 
[2] presented a paper at the International Symposium on Finite 
Element Methods in Flow Problems in January of 1974. They 
presented a method of solving problems related to floating 
bodies in water waves under unsteady-linear or stationary- 


nonlinear (second-order) conditions. 
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Eee OBJECTIVES OF PRESENT RESEARCH 

The experimental studies of Paulling and Porter [29], 
Paulling and Richardson [28], and Vugts [40], previously 
mentioned, indicated good agreement with the developed linear 
theory if the amplitude of motion was small. Potash [31] 
observed that small variations in amplitude caused large 
Giscrepancies between Vugt's experimental results and those 
of linear theory. Potash [31] has shown that second-order 
effects become significant at higher frequencies and are more 
significant when roll and/or sway modes are present. However, 
he did not examine the effect of finite depth on second-order 
eeluctions. 

, The present work has two primary goals. The first is to 
determine if the finite element method can be successfully 
applied to finding solutions to second-order boundary value 
problems of the type previously described. And, if success 
is achieved in the first objective, the second major goal is 
to study the significance of nonlinear effects in heave in 


water of finite depth. 
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f{1. FORMULATION OF THE MATHEMATICAL MODEL 


Wne pnysical problem is to determine the force required 
Mmemoustain steady, vertical harmonic motion of a floating 
Syeandric ship hull form. The ship will be considered 
eeerllating in an “ideal" fluid of infinite horizontal extent 


and arbitrary constant depth. 


Free Surface 








Bottom 


a, 7. mM 77a 


FIGURE 1. Hull at Free Surface 


A. THE NONLINEAR BOUNDARY VALUE PROBLEM 

The domain of the boundary value problem is shown in 
meee. In the figure, b is the ship's half beam, h is the 
mean fluid depth, dis the ship's draft, and 2A is the total 
Submerged cross-sectional area of the ship. 

The fluid is assumed to be incompressible, inviscid, 
homogeneous, and without surface tension. It is well known 
Maem unger the assumption of irrotationallity, a velocity 


potential exists. Furthermore, the potential will be a 
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single-~valued function in any simply-connected region. The 


@emvention for the potential used herein is 


Vo = ul + vj, | (1) 


where u is the x=-component of fluid velocity, v is the 
y-component of fluid velocity, and ~ and 4 are the unit 
vectors in the positive right~handed cartesian coordinate 
system directions x and y respectively. The symbol V is the 
standard "gradient" used in vector calculus. 

The nature of the problem dictates that at a large 
Mer. zontal distance from the ship's hull (cross-hatched area . 
fig. 1) only an outgoing one-dimensional traveling wave will 
be present. Therefore, it is advantageous to consider the 
meaweed region shown in figure 2. In figure 2, w is the 
region semi-~width. The meanings of the other symbols shown 


Mei ieure 2 will become apparent in the following discussion. 
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Radiation 
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FIGURE 2. Reduced Region 
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Conservation of mass requires: 
Vo? = 0 , ie eee ee) 


Ris the fluid domain and ve is the Laplacian operator. 
Bere are two nonlinear free surface conditions, one dynamic 
eeemone Kinematic. The dynamic condition, from the Bernoulli 


SPomation, is 
+ oo] +PR+yeg = OO. (3 
y 0 


Equation 3 must hold on the moving free surface (p = 0). In 
the equation, g is the acceleration of gravity, p is the 
overpressure (referred to atmospheric), ep is the fluid 
Seijoety, and t is real time. The subscripts denote partial 
differentiation with respect to the indicated variable. The 
kinematic condition comes from a basic assumption from 
Continuum mechanics simply stated by Stoker: "A particle 
once on the free surface remains on it." This assumption 


Geanslates into the following expression 


Lin ka) ae 


a ; (4) 


lin Sota ONE s.e unc pOssible pure function of time, 
which may appear, is assumed to be identically zero without 
loss of generality (See Stoker [35], p. 10). 
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Again, equation 4 must be satisfied on the moving free 
surface. In equation 4, D/Dt denotes the co-moving derivative 
(sometimes called the substantial or total derivative) and the 


function n(x,t) is defined by the following expression 
anew) ee (5) 


Therefore, n is defined to be the wave height above the line 
y = 0 at any point x and any time t. 
Equation 4 can be rewritten using the definition of 6 iyo 


equation 1 and the definition of the co-moving derivative® 


= + 
ony oy Ne 


a 0 3 (6) 


on the moving free surface. 


On and "near" the bounding surface Si (radiation boundary ) 
VG mene x yet) , as keto . oe 


In equation 7, #* denotes the velocity potential of a 
traveling wave moving away from the ship. It must be mentioned 
here that inherent in the statement of equation 7 is the 
assumption that the surface S, is "far enough" away from the 
ship's hull. This point will be discussed in more detail 
later. 

imrnane Co the interface So between the fluid and the 


Sap 's hull, the kinematic condition is 


rhe co-moving derivative is defined as: 


Meydt = ul + vo, + 0, 


eT 





Voen=Ven , (8) 


The vector n Hoe chew unic GuEward normal vector (from the 
meemon KR) to the surface at the point of interest. The dot 
is the standard scalar product and V ROpreosecnibom une Vecuor 
wemocity Of the ship at any instant of time. 

Since the motion of the ship is assumed to be restricted 
to the vertical direction (heave), the bounding surface So 


becomes a plane of symmetry. The surface 8S represents the 


51% 
mememom Of the fluid region and is rigid. Therefore, both 


Surfaces have the same kinematic boundary condition 
Von=0., (9) 


memexomination Of equations 3 and 6 indicates the nonlineari- 
ties inherent in the boundary value problem. An additional 
problem, not as obvious, is the fact that the free surface 
and the hull are moving boundaries. In its present form the 
meoplem has not yielded solutions. Therefore, some approxi- 
mate theory is applied. 

In the development which follows, a perturbation analysis 
will be applied which when carried to the second-order will 
produce two linear boundary value problems from the one 
Peesentily posed. Further, the boundary conditions of both 
linear problems will be referred to fixed boundaries. In 
weoeetom, appropriate equations for determining force, and 


Wave profiles will be obtained. 


C 8 \ \ 





BeeeeEERTURBATION ANALYSIS 
First it is assumed that @ and n each possess a perturba- 


tion series representation as follows 


od = ao ap oe + O(e 7) ; GoD 


eo sa i Pore + O(e >) ; (ia) 


where e€ is a small non-dimensional parameter which is a 
measure of the size of the ship's motion and is defined as 


meeeyb ("a" is the amplitude of the ship's motion and "b" 


mae (2) 


is the ship's half-beam). Also in equation 10, and © 


are defined to be the first and second-order velocity 


(1) (2) 


potentials, respectively, and similarly for n and 7 


Note that as e70 there is no disturbance and therefore there 
Peeemmo zero-order functions. Substitution of equation 10 


into equation 2 yields 
eyeo it) 4 62y2562) 4 O63) = 0, (12) 


Equation l2 implies that 


ooo) = 0) ane Gi 


vee?) = 0, anR. (14) 


Pesiintilar result holds for equation 9. 





ifmeeume Linearized Free Surface Condition . 
The following development of the linearized free 
surface conditions follows in general that of Stoker [35]. 
The final linearized equations for the free surface boundary 


emo .On are obtained in the following steps. Equations 5, 


mimeand 1] are substituted into equation 3, noting that p = 0 
Memeeeeceiree Surface, Also, it is observed that git), ’, 
and al le vavomre py resenvcamey some tOllowing form of the 
Taylor series 
o (x,y,t) - ,(x,0,t) + n(x, t)%)(x,0,t) 
né 3 3 
at = ® yy %05t) i (0p) ee (15) 


Micer collecting coefficients of like powers of e€ in the 


merrreeing equation, the following set of equations results 


(1) ~ ad = 0 (16) 


Bn 


27 gh) 4 nok) 4 Sealy? + (al }7) = oF, 


gn 5 


on S3 cy) 


Applying a very similar procedure to equation 4 yields the 


following result to the second order 


the subscript v denotes partial differentiation with 
mopeceu FO xX, y, or t. 
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et ae 


- . ona Cine) 
(2) (2) eG) (l)me(l) 
Ne - - + ie o - 7 os = 0, on S. (19) 


Differentiating equation 16 with respect to time, and 


eliminating n, from equations 16 and 18 yields the well 


t 


known linearized free surface condition 


go?) + gS) = 


bt > onye=odo., (20) 


Mmimeqiati1on 20, because of the use of the Taylor series 
Expansions in mn, the functions are evaluated for y = 0, 
Mme r than On the actual free surface. 

femmeeowing the same procedure, as before, for equations 


ieond 19 yields 


(2) Cee le tl (1) (1) 
eo + © = = >, (® + NY ) 


y cie g ey 
CA lh), Gil) 
= 2(o, ae a e Put ) : Can) 


a) 


Meme Lhe nonlinear function of 6 on the right-hand side of 
Eemwation 2l. 
Eeeelne Ship Interface Condition 
PnCe whens SH peboOundary is also in "motion, a 
Similar procedure must be applied there to "fix" the boundary 


So. The development follows to some extent that of Potash 


[31]. Consider figure 3 in which the general ship hull form 


Sel 





ax 


Ship Hull,S, 


PUGURE Ss. “Muli “Coordinate System 


is described in two-dimensional carteSian coordinates in 


terms of the parameter s (are length) in the following way 
paeex(S) 95). y = y(s) . (22) 


The inertial coordinate system is the same as shown in 
figures 1 and 2. The parametrically represented body 
coordinates (equation 22) correspond to the inertial coordi- 
nates of the body when the ship is at rest (neutrally buoyant 
position). The motion of the ship in this reference frame 


can be defined in the following way 
x(s,t) = x(s) , C23) 


Were ey) ete y(t) , (24) 


oye 





where the function Yn, is defined to be 


t 


= fRetee (25) 


Yh 


Meme aguation 25, o is the circular frequency of the harmonic 
Meuron and i is the imaginary unit. 
The velocity vector for the ship referred to the 


inertial reference frame is 
V=e hu 7 (26) 


mae aoc denotes ordinary time differentiation. The unit 
outward normal vector (from the region) in terms of the 


parametric coordinates becomes 


ne-y'i+ x'9 , (27) 


Mnere the prime indicates differentiation with respect to 
e@eevariable s. Substituting equations 26 and 27 into 


equation 7 yields the ship-fluid interface boundary condition 


t 


Voen = ey), x (28) 


Now expanding 6 in a Taylor series ine (similar to free 


surface expansion) yields 
@[x(s,t),y(s,t),t] = 0(X,¥,t) + eyo (K,¥,t) + O(e*) . (29) 


Using equations 10 and 29, the normal derivative of @ on So 


may be expressed as 
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— _ re elt —,, (1) 
Voren = (~-y oo + xX ae ve 


+ (-prel?? + x92) _ Fry O88) + ary, ol Je? + o(e3) , 


XY 


aie So (30) 


Pees elvuting equation 30 into equation 28 and equating like 
powers of e€ yields the following set of equations valid on 


the fixed boundary So: 


volt), 


+ 
N 


YnxX' > (31) 


(1) 


(2) 
Vor’. a 


+ Ty Oh) (32) 


Sy 
M 


an En 
Equation 32 may also be written 
Fol2) 2 = -y, (is¥o")? yy , (33) 


It remains to consider the boundary Sy (radiation 
boundary) and the function #*. For consistency, #* must 


also be expanded in a perturbation series 
ox = coed) 4 c2—2 62) 4 O63) | (34) 


Potash [31] and others have shown that px (1) and px 62). for 

a region of infinite depth, each represent a simple harmonic 
traveling wave of appropriate frequency and wave length. 
Zienkiewicz and Newton [43] have provided the appropriate 
homogeneous boundary condition for the surface Sy to Success-— 


fully pass a wave form without reflection. The equations are 
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Gox(l) 2 2 oxi) ; (35) 


all 
(Oe) op oe ae (36) 
Co 1 
In the above equations c, and C. (the wave celerities) are 


femened by the following implicit relationships 


= oh 
Cc, = (G:7 40m Man ve-hele! — ‘ 3) 
C5 = (g/20) tanh (£29 (38) 
O 


Seeeeeeruurbation Formulas for Pressure, 
Force and Wave Amplitude 


The determination of dynamic loads on the ship 
meauares the following equations which follow from the pertur- 


bation series assumed for @® and n 


Mey ,t) = p69) (x,y) v ep!) (x,y,t) + e*p*) (x,y ,t) t O(e°) : 
(39) 
Mee er 2) (4) + c@p'2)(t) + Ofe3) . (40) 


In equation 40, F(t) is defined as the total force per unit 
length required to sustain simple harmonic vertical trans- 
lation (see figure 3). 

Pollowing a similar procedure as before, evaluating 
the Bernoulli equation (equation 3) on the ship's hull using 


equations 23, 24, and 29 yields 
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: = p[x(s) ,¥(s) ,t] = gy(s) + eLgy,(t) + oY) (X96) 


+ "Ly ote + hee +i (901) 34 + (002) }*)) Soee) Chay 


By comparing equation 39 with equation 41 we find 


p = —- ogy(s) , (42) 

pi?) SOE reat o( 1); g (43) 
C2)ans (a) 1: Gin Gi (2) 

p == Dit eg = Bn. }= + ie }“) + o) ] 


(44) 


At this point observe that the second-order pressure is made 
memot contributions from both the first and second-order 
potential solutions. 

Potash [31] and others have provided relations for 
the force per unit length acting on the cylinder. However, 
memimis Lreatment, if 1s desired to obtain relations for the 
exciting force as previously defined. The motive is that the 
theoretical result obtained should correspond to that 
measured by an appropriate physical experiment. The relation 
between the exciting force F, harmonic displacement Yn and 


myarodynamic pressure p, is 


R 
F(t) = Ep2ay, - f pn: ds , Case 
-£ 


36 





where 2 is one-half of the total symmetric arc-length of the 
ship's hull (see figure 3). It should be noted here that 
Potash [31] has examined the question of the time-dependent 
total arc-length of the actual wetted hull and its effect on 
the limits of integration in equation 45. However, it will 
@emowaved Without proof that, provided the sides of the ship 
[emeeervaical at y = 0, the second-order effect 1s zero. The 
only hull form examined to the second order satisfies this 
stated requirement. Therefore, equation 45 is exact for the 
Case considered. 

iMmehc=carecullation of ~The force , the constant 


(0) 


mammeostatic lifting force which comes from p will not be 
memodered Since it is not in the definition of F. Accordingly, 
the following expression for F(t) results using equations 27, 
39, 43, 44, ana 45 
E % (1)\= 
F(t) = e[pe2Ay, - Jf -plgy, + o:°°)x' ds] 
—X 
+ e°[- Pe — Leah Ee “f cae aa < o 6?) yx ds J 
P Yh ay, 2 Ng y 1 ; 
ze (46) 


Comparing equation 46 with equation 40 yields the following 
set of equations 


ig 


P'Y) (t) = p2ai, + 0S Cay, + EET as, (47) 
L 
(e) p (1) rd le) 2 2 (2) 4=, 
fee (t) = 0 2 Lyn ty a ea fo. }) =e ix’ ds 


(48) 





The wave amplitude may be written by inspection of 


equations 16 and 17. 


a) Pg Oe, 
TN ee 2) 
m2 ¢x,6) = = 2606?) + nel) 4 Ee¢olt3? + 600132); 
(50) 


The boundary value problems for gf) and 962) may now 
Mmemeormally Stated. However, at this point, a standard 
separation of variables technique is applied to eliminate 
time from the solution (since the solution is assumed 
periodic). The appropriate definitions are 
lot, 


6) (x yt) = Pet) leans e : (51) 


Mitte) aRelo Gy) 6) (52) 


Note the newly defined spatial potential functions g 62) and 
g 2) are in general complex. Also because of the definitions 
in equations 51 and 52 it is necessary to replace ys, by the 


eomplex form 


y= be” (53) 
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C. THE BOUNDARY VALUE PROBLEM For 9‘? 
(1) 


The boundary value for @ PS coriaiiyestaved. Usine 


Seperrons 9, 13, 20, 25, 31, and 35. 


vot) = 0, ink , (54) 
gg 1) - 041) = (py eles (55) 
y 3 
Fo 6) 2 = 10m00' . on Sy (56) 
Be sey (1) 
Vo st). Hae » on Sy (oi) 
sail 
Pe 0 | on S, and S, . (58) 


D. THE BOUNDARY VALUE PROBLEM FoR 9 ‘?? 
(2) 


Similarly, the boundary value problem for $9 may be 


Meemally stated using equations 9, 14, 21, 25, 32, and 36. 


Poo? = 0, ink, (59) 
Ms ie2,(2) . 209°" -_42,(2) 4 (1) 
ey 26 y i 
so (60232 + £9192) | on ss, (60) 
=e) =, (1) 
by" 4 bi" ¢ 
$62) 2 = age = — yy » on Bs (61) 


the factor % on the right-hand side of equations 60 and 
61 is due to complex algebra involved and is explained in 
Appendix A. 
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: (2) 
Hy 62). aac aor Sy (62) 
Toned = 0, on S, and 5. (63) 


Note the striking similarities between the first and second- 
Sameer problems. Note also, that the nonlinearities in the 


Original problem have been approximated by the nonlinear 


(1) 


ferms in 


oL) 


in equation 60, and the nonlinear terms in he 


and $¢ imme diuacaom © 1 (Vy, does not explicitly appear). 


The additional nonlinear amplitude dependence comes from ae 
in the original perturbation expansions (equations 10 and 


iene The structure of the problem requires the solution of 


a1) Rey 


first in order to define the solution for 


The conversion: to complex algebra allows another observa- 


ee) 


tion to be made. Since only the real part of ¢ is 
meaningful, Appendix A demonstrates the required complex 


Manipulation. Close examination of that result indicates 


G2) 


that the nonhomogeneous -boundary conditions on 9 contain 


g 62) 


time-independent terms. This implies has a time- 
dependent and a time-independent part. However, since the 
M@aysical quantities of interest in the solution all involve 


(2) the time-independent solution of be) 


derivatives of $@¢ 
will not be pursued further.°> The above comments are also 


meue ior the second-order functions of pressure, force, and 


Lee [18] comments on the time-independent solution in 
terms of mass transport. 


HO 





wave amplitude. However, the time-independent quantities 
are of interest and have physical meaning as will be seen 
later. 


(2) 


The boundary value problem for $6 is valid only in 
\ 

water of infinite depth. The following discussion will 

explain the problem and develop a boundary value problem 


cay 


Meech will allow for the solution of 6 in finite depth. 


E. THE SECOND~ORDER BOUNDARY VALUE PROBLEM 

Wel FINITE DEPTH 

Potash [31] has shown that the asymptotic form of the 
velocity potential @ as x>* becomes (in water of infinite 


depth) 


i Lote (ky x+y i[20t~(ukyxty‘*?)] | 


(»,0,t) = Re{eG,e + e°G,e 


(64) 


woere Ky denotes the wave number of the first order wave 

A 
(ky = o/c) and Gy and Go are coefficients which depend on 
frequency o. However in the more general case of finite 


Pepun, an additional term due to the Siecieas 0” second-order 


potential appears as follows 


F (1) (2) 
d(>,0,t) = Re{eH,e bot (k1xty )] +e°[H,e 1 Katy ) 
O 
al 0, oe se or2i(kyxty <2), acist, ; (65) 
Bc, sinh (k,h) 


Ome Stokes potential to the second order is well known 
wo be zero in infinite depth. 


Wy 





fmarequation 65, k, is defined as the second-order wave number 
(k, = 20/c 4). Hy and H., are constants similar to G) and Gy. 
Close examination of equation 65 reveals that the radiation 
Bemndary condition (equation 62) would fail since the 

"Stokes wave" which is present is traveling at the same 
celerity as the first-order solution. This means that the 
radiation boundary would "reflect" the Stokes! wave since 

the boundary condition defined by equation 62 will only pass 
a wave traveling at a celerity of c,, defined by equation gS. 


A successful resolution of the above described difficulty 


requires the formulation of a different boundary value 





problem. We define the complex velocity potential y 62) as 
2 
(2) _ azo TH coshleky (yt) ag (ay xty)?) P 
ee 2 tl, (65a) 
Bc, sinh* (k,h) 


which corresponds to the third term in equation 65. Then, 


following previous conventions, the complex velocity potential 


72) 


is defined as 


1. The Second-Order Boundary Value 
Problem for 6 - 


The definition of G2) in equation 66 allows the 


formal statement of the following boundary value problem, 


(2) 


which follows from equation 66, the definition of yp and 


(2). 


the boundary value problem for $@ 


he 





72 562) = 0 


» ink , (67) 
ep 2) — yo2Gl2) 2g 99? 65241) 4 (2); 
y 26 y Ee yy 


os tof {g<)}° +f (o61) 32) - eve? + ug 2y 62) , on S 


Si 
(68) 
Ca) (1) 
, (be _ b6 
Ho 62) sal (— + yi??) = xt (Se a bee?) » On So» 
(69) 
; ee ( 20) 
vo?) 3 = = , on Se (70) 
ieeen = - Ws -n , on a Gal 
Woe) nl =O On S. (72) 


Equation 70 will now properly behave because the asymptotic 


=( 2) 


form of uSing equations 65 and 66 is 
; G25) : 
—(2 Aa) eG —~2i(k5xt 2 
GP) (oy) ertot - Hy e pase) sees (73) 
Note the celerity of the complex wave form gi ecist is Co 


mae tne radiation boundary condition will now work. 


FF, FORMULAS FOR PRESSURE, FORCE, AND WAVE AMPLITUDE 


The development which follows requires additional defini- 


tions, necessary because of the separation of variables scheme 


previously introduced. They are: 
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(1) tot, 


ps2) (x. y,t) = Re{p ‘ p62) = Re (pte? erto®} | aly 


Meee) = Rete 7e2O*) | eA) (4) = Re{ri2 ctor, (75) 


nit) (x,t) - Re {7 $2) (xyet%*} ; n 62) = Rel CU Gee. (76) 


(1) =-(2) =(1) =(2) 


The functions p cals a Ri) are complex. This is 
also true for the constants fl) and phe) 
Using the above equations and equations 43, 44, 47, 48, 


49, and 50 the equations for pressure, force, and wave 


Memene for both first and second-order solutions follow, 


pet) = -p[pg + tog] , (77) 
reer. ) 
jog s~?b 
pee? = pp F + S92? + (909?) + atop], (78) 
L 
Mee —>opac? + p S (ep + tog‘ ))x as , (79) 
-2 
Ci) 
J = alent b a 
2) -p Jf (—3— + 5 ae a a (962) }°3 + 2iag\*) yx ds, 
-2 
(80) 
Ci) 
7 (i) a —, ; (81) 
/ Gee (ier 
. —(1),(1) {¢ } {¢ } 
=(2) = - = (2tgg{*) + ton” ty + Sp + 1) 


2 
(82) 


Again, equations 78, 80, and 82 include the special considera- 
tion of the complex algebra involved which is explained in 


Appendix A. 
Wy 





eee ELEMENT DISCRETIZATION 


the Finite Element Method and its application to scalar 
mend problems is described in various texts, including one 
by Zienkiewicz [45]. The essential ideas involved follow. 
m@emrepion of figure 2 is replaced by a grid or mesh of 


finite elements as shown in figure 4. 





FIGURE 4. Finite Element Mesh 


Mem it 1S assumed that the field function sought can be 


represented by the following equation 
o(x,y) =N o> , inR. (83) 


The total number of nodes in the region will be defined to 
be m. Then N is a 1 xm row vector of shape (interpolating) 
functions and ¢ is an mx1 column vector of nodal values of 
the complex potential ¢ (in this problem). Equation 83 
implies that knowledge of the vector 9 defines ee. Vo) 


Everywhere in R. 
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inere are several choices of shape functions Ny which may 
be used in equation 83. This work selected the cubic shape 
functions described by Zienkiewicz [45] as the "serendipity 
family". The elements used are isoparametric. A brief 
description of the l2-noded isoparametric elements is given 
in Appendix B for the unfamiliar reader. The key advantage 
of the isoparametric elenene 2 lessin they tacy that regions 
with curved boundaries may be readily represented. 

The development that follows defines the set of linear 
equations which determine the vector $ of equation 83. 

One general approach to discretization of several types 
of field problems in finite element analysis is to define a 
functional whose integral over the domain of the problem is 
meme minimized. The Calculus of Variations shows that if the 
Euler equation of the integral to be minimized is the governing 
equation (i.e., Laplaee's Equation) then the two problems~are 
Peiivalent. 

However Zienkiewicz and Newton [43] have shown the 
discretization may be readily accomplished by applying the 
Galerkin process. The weight functions are chosen to be the 


nodal shape functions. Mathematically stated 


on A OL ETRE ta Oe (84) 


This statement assumes the same boundary conditions are 
moed. 
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Meee superscript T denotes transposition and the function 6 
mamconsidereda to be any of the complex potential functions 
Maese SOlution is sought. Application of the divergence 
theorem transforms equation 84 into 


sVo-n as, (85) 


T 7 
S €NO No] fo aR = fsN 
oy 


where S denotes the boundary of the region R. 
A. THE DISCRETIZED FIRST-~ORDER PROBLEM 
Application of equations 84 and 85 to the first-order 


boundary value problem (463) ) stated in Chapter II yields 


mae TOollowing matrix equation 
a, +22 p +H) og’)? = tone)? | (86) 


The definitions of the matrices in equation 86 follow. The 


meee 1 1S mxm, real, symmetric, and defined as 


a T OT : 
H : [Ny Nid N,| aR . Co) 
Pay 


The matrix H comes from the integral on the left-hand side 
81 equation 85. The matrix D is mxm, real, symmetric, and 


defined by 


T (88) 


Vo 
il 
on 
eas 
Ue, 
(er 
op) 





Sees CONtCribuced by the boundary integral on the right-hand 


Side of equation 85 and represents the discretization of the 


radiation boundary condition (equation 57). 


The matrix Q@. is real, mxm, symmetric and given by 
i WOW else C32) 
3 an 


Qo parallels Din its origin and represents the contribution 


Sd 
~ 


from the homogeneous free surface boundary condition (equation 


Bp). 
co) 


The vector r 1s real and defined by 


pit) = f Ine seus shee (90) 
e G 


2 


The above vector is the discretized ship-fluid interface 
kinematic boundary condition (equation 56). The surfaces 


Si. S35 Sys and S65 have homogeneous boundary conditions and 
therefore make no contribution to bl). fp iol ows chat ae? 


is non-zero only at nodes along the ship's hull (surface So). 


B. THE DISCRETIZED SECOND-ORDER PROBLEM 


se) 


The problem for $6 only will be presented since it 


merivains the solution for g 62) Pic nceLinitu. tne boundary 
(2) 


value problem for 6 yields a linear system of complex 


equations similar to equation 85. 


(-40*Q. + Ee p 4 H) g 62) = (2) (91) 


oe. 





2) 


The right-hand side vector r in equation 80 receives 


Senera butions from the surfaces Sy > S5 and D3: Reeoradinieiy 


a further decomposition is necessary: 


2) en 2) 2) (92) 


x Sil u2 ae 


cz) 


mae Vectors rs feemec mile G sel Ow.  jubrst to be COnsistenc, 


(2) 


mae function y 62) must be represented in the same way as 6 


G2) 


Then ry becomes, from equations 71 and 85, 


ie =f NN (2) as. (94) 
se g ies 
a 


(2) 


- are determined by partial differen- 


C2) 


The nodal values of w 


meauvon Of the function w and then representing the 


(2) 


x in similar manner as equation 93. 


mone TION wy 





The vector eee? is defined by equations 69 and 85 
lading (Dp 
$2) = 5 ntny (=X — + yl) ya5 - s wnxt(=¥¥— + yl? Jas 
x g ~~ x q = = ~y 
e (95) 


Sima barly Se follows from equations 68 and 85 


(2) 1 T . 4g 40 (2) 2 (2) 
Meee ee EF OY By toys 


3 (96) 
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The ee component of a is 


Ga Gly) 


I a (97) 
where p61) and oy? slides SW IGg\= ae components of gt) and bo, 
mespectively. Similarly 
Rey) 2 (1) 
By = 70° Oy yy » and (98) 
Pein (99) 


as) 


The nodal values of the various partial derivatives of $¢ 
are determined by differentiating equation 83. The method 
Seca lculation is discussed in the next section. The reader, 
PMaamiliar with isoparametric elements, should refer to 
Appendix B prior to reading the next section. 


te Caleulataeon of Nodal Values of the 


Field Derivatives of $¢ - 


Ci g 61) gf) 


required in 
deer yo CY q 


The vectors $6 
equations 95 and 96 are determined by applying standard 


Bemcepts of finite element analysis. The appropriate relation 


th 


mer the i node of a given element is 


(l)e e 
Os Xx N 

CIDE ~ e ~ 
Ory Nn 
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Jv is the 2x2, real Jacobian matrix of the transformation. 


The vector N© is a 1x12 row vector of element level shape 
functions. The superscript "e" denotes all quantities in 
Brackets are at the element level. In equation 100, the 
Ewements of ee and N° are evaluated at node i. 


At 


Application of equation 100 to all twelve nodes of a 


given element determines the vectors va and _* ° 
Then, replacing g lle by at in equation 100 yields the 
components of the vectors grve and due Mathematically 
Stvacved 

ae NE 

= gre poe (101) 
gfe ye 
15YV | ~ 
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IV. THE COMPUTER PROGRAM AND NUMERICAL SOLUTION 


Pee MAIN FEATURES OF THE COMPUTER PROGRAM 
iice.Ormabion and Solution of the linear system of equations 


defined by equation 86 Cie (2) 


system) or equation 91 (6 
system) was accomplished on an IBM 360/67 computer at the Naval 
Postgraduate School. The computer has a maximum word length 
of sixteen bytes. 

A comparison of equations 86 and 91 demonstrates the close 
Similarity of the coefficient matrices for the first and 


second-order solutions. Note that only the coefficients of 


25 and D change. An additional observation follows from 


we 
ed 


examination of equations 87, 88, and 89 defining H, D, and 

Qo respectively. Note that all three matrices are purely a 
function of mesh geometry. Therefore, once a mesh is chosen 
memeeetine the region of interest, the three matrices previously 
mentioned need only be calculated once regardless of the forced 
moui1on . 

A well known property of the finite element method is that 
ie wsually produces a banded, symmetric, coefficient matrix. 
That is the case here and this has obvious computational 
mimanvases. However, in this problem, the final system of 


equations is a function of the frequency of motion as for 


example in the first-order problem 


ow 


Sz 





(1) 


The matrix K is an mxm complex matrix defined by equation 


poco be 


K{*) = -0%Q, + 


we 


Bl se ie cles) 


Pac hough Bo? Pocomelex ana frequency dependent, it is 
possible to perform most of the triangular decomposition of 
ee? (by Gauss elimination) only once, using real arithmetic, 
for a given problem geometry. This economization was first 
employed by Chenault [4] and is described in more detail in 
Peepemaix C. In addition, this author employed a vector 

Cae) 


Storage of the upper triangular portion of K 


~~ 


reduced the storage requirements of the program. The vector 


we wie Grea talk, 


storage scheme made it possible to store only the portions of 
ai which are non-zero. 

The net advantage of employing Srenecle nee scheme and -the 
compact vector storage was that the system of 733 complex-. 
equations which is generated by the mesh shown in figure 4 
could be solved for both the first and second-order problem 
twenty times in twenty minutes. Further, only ten percent of 


the total matrix was actually stored and the computer solution 


was accomplished in core. 


Bone Chenault scheme described in Appendix C was developed 
by Professor R. E. Newton at the Naval Postgraduate School. 
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B. NUMERICAL SOLUTIONS 
The solution of the first-order problem for a given mesh 
geometry is required before the second-order solution can be 
Memeined, Purther, the first-order solution must be highly 
accurate if a good second-order solution is to be expected. 
This fact is made more obvious by re-examining equations 68 
and 69 in part II. 
meee irsy Order Numerical Solutions 
several first-order problems were solved and closely 
examined for accuracy. The numerical solutions are presented 
mpart V. 
a. Mesh Requirements 
One of the primary considerations in choosing a 
Peeper mesh concerns the positioning of the radiation 
boundary (S)). It may be recalled that the surface Sh must 
be far enough away from the ship hull (location of disturbance) 
so that the local effect of the disturbance has decayed and 
only an outgoing wave is present. Bai [3] has addressed this 
question analytically. He has shown that the eigenfunctions 
associated with the non-propagating velocity field decay 


exponentially with x. Bai recommends the criterion 


"agent D1 (104) 


where ais a monotonically increasing function of frequency. 
The range of ao is from 1.5 (low frequency) to 3. The values 
of a were determined based on the decay rate of the most 


persistent eigenfunction and an attenuation factor of 0.01. 


5 | 





This investigation has found the above criterion somewhat 
Bemservative. A value of a = 2/3 was found to be adequate 
meemune Velocity potential to numerically approach its 
asymptotic limit for all frequencies studied for the semi- 
Meeecular hull (h/pb = 6). Similarly, a value of a = 9/10 was 
found to be adequate for the bulb hull form described in 
part V (h/b = 10) for all frequencies studied. The values 
of a specified above were determined for the deepest regions 
studied for each hull form. The region widths used in this 
case (w/b = 6, w/b = 10) were found to be more than adequate 
for shallower depths. 

Another requirement on w which tends to restrict 
its maximum possible value comes from the practical limits on 
fare ineness and the resulting computer storage requirements. 
The mesh must be sufficiently fine to represent properly a 
weaver surtace wave. The finer the mesh, the largertethe _ 
computer storage requirements and the longer the computer 
wemutrton time. Therefore, this latter restriction tends to 
oppose the former one. 

Bai [3] used quadratic elements and recommended 
at least 10 surface nodes per wave length. Visser and Van 
der Wilt [39] also used quadratic elements and recommended 
that an element could not span more than one quarter of a 
wave length. This investigation examined this question in 
Seca. As previously mentioned, cubic isoparametric elements 
were used. It was found that if tne gravity wave did not have 


to be propagated more than two wave lengths, the elements 
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could be allowed to span six-tenths of a wave length. 
However, as in the case Visser and Van der Wilt [39] studied, 
the successful propagation of a gravity wave over many wave 
lengths requires a finer mesh spacing. The one quarter wave 
length Visser and Van der Wilt recommend was found to be a 
Meactical upper limit. 

A parallel problem to positioning the radiation 
boundary concerns the location of the bottom (surface S.) 
to simulate infinite depth. Comparison of the finite element 
memuetons obtained in this study with analytic solutions 
based on infinite depth led to the conclusion that a keel to 
pottom clearance (h — da) of 5b will adequately simulate 
Mmanice depth. However, satisfying the classical hydro- 
dynamic requirement that the depth must be at least one-half 
of a wave length becomes impractical at low frequencies. 
Ignoring this requirement did not seem to affect solution 
meecuracy. 

b. Mesh Data Preparation 

The meshes for a given boundary geometry were 
prepared using a very versatile mesh generator computer 
program developed by Adamek [1] at the Naval Postgraduate 
School. This tool greatly expedited mesh preparation. For 
example, the entire data for the mesh shown in figure 4 
(733 nodes, 132 elements) could be generated in twenty 
seconds. A modification was used to supply accurate values 


of the nodal coordinates along the ship-fluid interface. 
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ec. Convergence of First-Order Solutions 
Convergence of the first-order solutions could 
not be determined by standard numerical techniques (i.e., 
mesh refinement) due to computer size limitations. However, 
the excellent agreement between the first-order solutions 
obtained in this study and both theory and experiment left no 
femby that the solutions were quite good. 
ew ocecond—Order Numerical Solutions 

The successful Sar eaiO of the first-order problem 
made study of the second-order problem possible. 

a. Mesh Requirements 

All of the requirements on mesh geometry previously 

mentioned for the first-order solutions of course apply to the 
eeeond—-Order solution. However, the second-order problem 
demands even more stringent requirements. One. reason for this 
follows from the fact that the wave lengths of the propagating 
portion of the second-order velocity potential are .much 
shorter (one-fourth of the first-order wave lengths in infinite 
depth solutions). Therefore, since it is impractical to use 
Pei terent mesh for ¢€ach solution, the mesh must have a free 
surface element spacing approximately one-fourth of that 


(1) 


required for the solution of $6 for a given frequency range 
to be studied. The mesh shown in figure 4 was used to solve 
the infinite depth second-order problem. It contains 132 


elements and 733 nodes. 
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pee Ine seeeno—CUrder Boundary Condition on 53 
Additional constraints on mesh geometry were 
found to be required. These constraints arise from the 
requirements to calculate second-order partial derivatives of 
g(t) which appear in equations 95 and 98. It was essential 
Mmemaave rectangular elements with equal node spacing along 
the free surface for success in representing the non- 
homogeneous boundary condition given by equation 96, 


theme ner alscusslOn of This DPOint will be 


aided by the definitions of some complex functions. Define 


ey (1) 
Q(x,0) = S28 — [0792 + gol] — solo a® + (6077977 , 
(105) 
feo) = gv”? ~ ygt ys?) , (106) 


In terms of these definitions the right-hand side of equation 
68 is 
P(x,0) = Q(x,0) - R(x,0) . COT) 


The reader should recall that the successful solution of the 


2) 


boundary value problem for 9 using the method of finite 
elements required the definition of a new boundary value 
problem for the function g2) The reason was that in the 
finite depth problem a second-order Stokes! wave is produced 
moecgin celerity Cc, as well as a second-order wave with celerity 
Co (contained in gle), The presence of the two waves with 
aGifferent celerities meant the radiation boundary condition 
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would fail. The problem was mathematically eliminated by 
Subtracting the well-known analytic solution for the Stokes 
wave portion of the second-order solution of g 62) One Very: 
boundary of the region R. The mathematical steps are straight- 
forward as shown in Chapter II. The asymptotic value of the 
Mmmet ion P(x,0) is zero as x gets large. But, the function 


R(x,0) which represents the classic Stokes second-order wave 


free surface boundary condition has the form 


5 ; (1) 

2 - + 

R(x,0) = a [eS ey) | oo EES ay, =. AO.) 
mer ita kd 

everywhere on 53: Micrelere. ti sehe numerical representation 

6f the function P(x,0) is to approach zero near the radiation 


boundary; the numerical representation of the function Q(x,0) 


must approach R(x,0) or 
OeuwO sak ( x0) asm (109) 


(1) must be known to 


(2). 


Note that the amplitude and phase of 6 

completely define R(x,0) and therefore jp 
The above discussion points out the necessity of 

a very accurate numerical determination of the function Q(x,0). 


—(2) 


The second-order solutions for $@ were closely examined near 
the radiation boundary and it was found that on the average 
the function Q(x,0) was within one-percent of the analytic 


fmemetion R(x,0). This accuracy was achieved even though 


formation of Q(x,0) requires calculation of second derivatives 
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Ot gb) using element shape functions which guarantee only 
C° continuity across inter-element boundaries. 


ec. The Ship-Fluid Interface Second-Order 
Poundary Conds ton 


The accurate representation of the ship-fluid 
interface (S5) boundary condition for the second-order 
boundary value problem (equation 69) is of course essential. 
Potash [31] encountered difficulties in making a numerical 
calculation on this boundary because the potential singulari- 
ties which define his solution lie along So. This results 
in the required functions pi? and are being only piecewise 
continuous. This problem required an additional approximation 
eo be made. 

The method used herein also has the problem that 
—- = bee and . are not continuous across element 
boundaries. However, no singularities are present. Extensive 
investigations in this work have shown that it is essential 
that the elements along the hull be as square as possible, 
and further that the node spacing on the elements along the 
fee should be uniform. It was found that if the hull 
elements were distorted (not square or rectangular) the 
numerical representation of the functions be and boo was 
highly unstable. However, the mesh requirements for 
momerically calculating pg? and oo? alone Tie vault were 
much less stringent. 


The adequacy of the hull elements used in this 


investigation (see figure 4) was tested by defining a known 
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mevential function on the elements along So 


le sy al sey (110) 


and then applying the computer program's numerical scheme 


to calculate the functions yy and boy: Observe that 


Pay aot. wae Cla» 
and 

t — 

boy ilies (112 


It was found that the cubic, isoparametric elements along So 
merece d an average error of about one percent in calculating 
fee tunctions byy and a Piecw iia xinuUMscmrer ODServed was 
five percent at isolated nodes. 

ad. Convergence of Second-Order Solutions 

Again, computer capacity limitations prevented 

the application of systematic mesh refinement to determine 
Bemverpence of the second-order solutions. However, it was 
possible to refine the elements next to the ship's hull and 
no significant change in the solutions was observed (the 
elements along the free surface were performing properly). 
mere: , If 1S Shown in Chapter V that second-order results 
for infinite depth are in satisfactory agreement with those 
of other researchers. No comparison is available for finite 
a ieiecoMmMbbons bub a Smooth transition is observed as the 


depth becomes effectively infinite (See Chapter V). 
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V. NUMERICAL RESULTS 


A. DEFINITIONS 

The definitions which follow are useful for presentation 
eeeresults. 

1. Added Mass and Damping Coefficients 


Mae of the first-order 


The complex amplitude F 
resultant hydrodynamic force acting vertically on the ship's 
mid is 


Q 
p(l) = 5(l) ss ds . (113) 
y —~2 


(ly 


fee real part of By is in phase with ship's acceleration 


and the ship's added mass is defined by 


Re{Fy} 


i) 5 Onno ae (114) 
o-a 


where the ship's acceleration amplitude is given by oa, 
The added mass coefficient is defined to be the ratio of 


added mass to ship's displaced mass 


Mm, Re {Fy} 
oe = SOA = oe : CLS, 
P ooAac 


jhe imaginary part of mt tSeiinomase Wich the 
9 


velocity of the ship and is the "non-conservative 


"The potential theory is, of course, conservative. However, 
Bae waves Which propagate to “infinity" represent work done by 
the ship which is effectively lost. 
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Semponent. Accordingly, the damping coefficient is defined 
to be 


~ImtF*)} 


Cc, = ———— (Lane 
é 2phac* 


Cn and Ca See Eeaumency dependent Tor a given hull form and 
region depth. 
A dimensionless frequency parameter is defined as 


follows 
2 


§ = c= (117) 
The parameter 6 is commonly used by other researchers and is 
used herein for ease of comparison of solutions. The 
parameter 6 is equivalent to eb nih oe o haga ols Oh oX- Wao (- Soh) ol @7, NEE a ol 
wave length). 
The added mass and damping coefficients previously 
defined are for heave. The same coefficients for sway motion 


Bee mear ty defined. The only difference in the definitions 


(1) 


is pil) is replaced by the horizontal force amplitude ee 


J 
defined by 


QR 
pi }) So % gp Jan be. (118) 
—% 
where 561) is the hydrodynamic pressure due to the swaying 
motion. There is a coupling between sway and rolling motion 
for a general ship cross-section. However, the only hull 


form studied in the sway mode was semi-circular in 
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Sross—section and the coupling force for this case is zero. 
Therefore, the appropriate coefficients typically used are 
memecetined herein. 

One further comment is warranted here. The second- 
order time-dependent hydrodynamic forces are harmonic with 
frequency 20. Therefore the definition of second-order added 
mass and damping coefficients is meaningless. 

e. Dimensionless Force Coefficients 


(1) and p62) 


The force coefficients f are also 


presented in dimensionless form. The non-dimensionalizing 
factor for the various force amplitudes is 1/pep’. Table Il 
presents the correspondence between the dimensional and non- 
dimensional force quantities. The second-order quantities 
tr?) ‘and ee which appear in table I are amplitudes of the 
time-independent (static) and time-dependent (dynamic) 

(2) (2) 


respectively. The existence of f WS. 


mermcions of f i 


Beigaiyeverified by examining equations 75 and 80 in light of 


Boe Giscussion in Appendix A. 


TABLE I 
DIMENSIONLESS FORCE COBFFICIENTS 


Force Amplitude eit) oe 
Dimensionless pCi) a nee 


Force Amplitude 


64 





3. Wave Amplitude at Infinity 
Another suitable measure of the damping of ship's 

motion (to the first-order) is the wave amplitude at 
Bemeimivy™ (i.e., near the radiation boundary). Further, 
experimentally it has been found to be an easier quantity 
to measure (particularly for high frequency motion since the 
damping force is very small then). Accordingly, the dimen- 
Sionless first-order wave amplitude at infinity is defined 


by 
Can) 


B. HULL FORMS 

Other investigators have studied a variety of hull forms 
Both by experiment and analytically. The semi-circular hull 
was chosen because of the wealth of data available for 
momearssen Of both first and second-order solutions. The 
bulb-shaped hull studied by Paulling and Porter [29] was 
also examined for first-order solutions. The shape of the 
Pulb hull in comparison to the semi-circular hull is shown 
ma figure 5. 

The coordinates of the bulb hull studied are defined by 


the following mapping 


e 
Re aye c+ ese (120) 
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The 


Centerline 






— Waterline 


Bulb, d/b=5 


Semi-Circle, d/b=l 


HAGURES eeu: Forms 


eomolex mumber c represents the points on a circle in 


the t-plane. The choice of radius in the z-plane defines 
the scale of the hull cross-section. The values of the e, 
are given in table II. 
AB Bele 
CONSTANTS FOR HULL FORMS 

Sioae ey C3 Te A/bda 
Semi-Circular 0.0 6.20 0.0 0.7854 
Bulb -~0.7132 +10) SO 2eralé 0.0605 0.6950 
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fee PE RST—ORDER NUMERICAL SOLUTIONS 
ie oemi—Circular Hull 


The semi-circular hull form has been studied by many 
meeearchers mentioned previously in part I. Specifically, 
eae Original solution by Ursell [38] for heaving motion has 
been confirmed analytically by Porter [30], Vugts [40], and 
@eners. Ursell's solution has also been experimentally 
confirmed by Paulling and Porter [29] and Vugts. The 
Sie eyone results of Ursell for Cr and Cy are shown by the 
Semr nucus curve in figure 6. The numerical results obtained 
by the finite element method (FEM) in this work are repre- 
sented by the circled points in figure 6. The FEM points 
shown were obtained using a mesh similar to that shown in 
figure 4, but using only 62 cubic elements and 359 nodes. 

iiomenvenoLton Of the Selutions for Cn and Ca TO 
finite depth is easily accomplished by simply redefining the 
mesh geometry (decreasing the vertical dimension). Kim [14] - 
and Yu and Ursell [42] reported analytic results. The 
results obtained by the FEM in this work for h/d = 4 agree 
meri those reported by Kim for Ci and nil? Wat isietd ee ieee 
percent over the range 0.25 < 6 < 5. The results of Kim are 
in disagreement with those obtained by Yu and Ursell for 
h/d = 2. The results obtained herein fall between the two. 
fees reports inability to verify Kim's result, but gives 


no details. Bai used the finite element method to obtain 


fees SOlUtiOns as mentioned in part l. 
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FIGURE 6. Added Mass and Damping Coefficients in 
Heave, Semi-Circular Hull, Infinite Depth 


| 


Figure 7 shows the effect of finite depth on Ci 
Sevained in this work. The figure shows that significant 
emanees dO not occur until very shallow water is encountered 
Maya < 2). Further, the effect of shallow water on C,, is 
more pronounced at very low and very high values of 6. Note 
the unusual variation with depth for 6 = 0.5. om decreases 


with decreasing depth (opposite the trend at higher values 
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2 ot eee 


FIGURE 7. Variation of Added Mass Coefficient with 
Depth, Semi-Circular Hull in Heave 
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of 6). Although it does not appear true in figure 7, oe 


(tor 6 = 0.5) has reached its asymptotic value at h/da = 6. 
Figure 8 shows the corresponding effect of finite depth on 


C The figure indicates that Ca is only sensitive to 


a 
Shallow water effects at very low values of 64. 


1.25 





FIGURE &. Variation of Damping Coefficient with 
Depth, Semi-Circular Hull in Heave 
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ihe semi-circular hull was also studied in sway 


motion to the first-order. The only changes in the boundary 


(1) 


value problem for @ defined in part II are in equation 56 


which defines the surface Sy: The changes are 
Gost) 2 = iobx' , on S.5 ; (125) 


for equation 56, and 


d = 0-5 om Ss ; (i222) 


momeeduation 53. The boundary condition on S. defined by 


) 


equation 58 remains unchanged. The ship's motion is defined 


in a manner analogous to equations 23 and 24 


mS.) x(s) + € x(t) : (225) 


VCS 5a) SSE (124) 


The function xX, (t) is defined the same as yy, (t) 


Oo 
= Re{be “} , (125) 


Tp 
and e€ is defined in the case of sway motion to be 


(126) 


02) 
i 
o|m 


The symbol "a" is the amplitude of the sway motion in this 
case, | 

Vugts [40] presents both experimental and analytic 
melutions for Ca and C, for the semi-circular hull in sway. 


d 
Figure 9 shows a comparison of Vugts! theoretical solution 


Heal 








FIGURE 9. Added Mass and Damping Coefficients in Sway, 
Semi-Circular Hull, Infinite Depth 
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with the FEM results obtained herein. The FEM calculated 
points are shown by the solid circles and squares in the 
figure. The solid curve (from Vugts' theory) was scaled 
from a very small graph and may be somewhat in error. 
meee BULD Hull 

Paulling and Porter [29] reported both analytic and 
Peoerimental results for the bulb hull form. They show 
excellent agreement by comparing first-order exciting force 


amplitude coefficients Fi) and wave amplitude coefficients 
Ad) 


wl) 


respectively. Figures 10 and 11 compare the values of 
and Ait) obtained by the FEM herein (circled points) 
and Paulling and Porter's theoretical calculations (solid 
line). The FEM mesh parameters were h/d = 2 and w/b = 25. 
This value of w/b was found later to be quite conservative 
(w/b = 10 would have been adequate). 

Figure 12 shows added mass and damping coefficients 
for the bulb hull vs h/d for various values of 6. The curves 
meemaesimilar trend to those of the semi-circular hull form. 


5(1) 


ihe corresponding curves of force and wave amplitude 


Ai)? are not presented. However, the maximum change in 
exciting force amplitude does not exceed five percent for 
hyd) = 1.2 as compared to infinite depth (h/d = 2). Similarly, 
the maximum change in abt) was only ten percent over the same 
range of h/d. 

Bee wave Generator Studies 


MaceanveleOln adn 19575 reported an analytic theory 


using a Green's function approach for solving the problems 


is 








—— Theory (29) 


FEM 
O h/d = 2, w/b= 25 


MPEGURE 10. Exciting Force Amplitude in Heave, 
Bulb Hull, Infinite Depth 


— Theory (29) 
0.4 FEM 
O h/d=2, w/b= 25 





FIGURE 11. Wave Amplitude for Heave, 
Budloe wue eee tind Ce Depth 
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h/d b/d 


FIGURE 12. Variation of Added Mass and Damping 
PEL Coefficients with Depth, Bulb Hull in Heave « 


Of pravity wave production by a moving partition and gravity 
wave reflection from a horizontal strip. The problem of 
gravity wave production by a moving partition was studied in 
this work to further help determine proper mesh characteris- 
mes Tor gravity wave propagation. 

The wave generation problem is readily placed in the 


framework of the class of problems being studied herein with 


1) 





10 


Soy minor modifications. All of the boundary conditions 


(1) 


Pesed for 6 in part Il remain unchanged except those on 


Sy and So + The surfaces S3 and So become a vertical line. 
The region therefore is rectangular. The left-hand boundary 


will be designated as Sy as Shown in figure 13. 


K Free ete 5 
“7 Wave Generator, 


h ¥ (flapper type) 
| / /-Bottom,S. 


FIGURE 13. Reduced Region for Wave Generator Problem 


“Ty “a 
/ 


Radiation 
Boundary,S, 


The boundary condition on Sj becomes 
oF = aly), on st, (127) 


The function q(y) describes the mode of displacement of the 
wave generator. The plunger type and flapper type (shown in 


figure 13) wave generators were studied and compared with 


o This observation is not true in the potential theory 
Semeutearity distribution approach. The source singularity 
must be replaced by doublets and the problem completely 
reformulated. 
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MacCamy's results. The function a, (y) for the plunger type 


generator is simply 


a, (y) Sus (128) 


where Ho is the maximum amplitude of the wave generator 
motion. Similarly, the function dp(y) for the flapper type 
Benerator is 

Cases) 


ey) = pe (iz? 





The comparison of MacCamy's results with those obtained 
herein is shown in figure 14 for the plunger. The circled 
points were obtained by the FEM and the solid line is MacCamy's 
Mmeoretical curve. Similarly, the comparison of MacCamy's 
result for the flapper type wave generator vs the FEM results 
~ obtained in this work iS shown in figure 15. The excellent 
fsrcoement allows two observations to be made. First, the FEM | 
approach to solution of problems of the type being presented 
is highly versatile in its application. Second, the radiation 
boundary placement was further verified. 
D. SECOND-ORDER NUMERICAL RESULTS 
1. semi-Circular Hull — Infinite Depth 
The first-order solutions just presented established 
Pomtvdence that second-order solutions could be attempted. 
The only hull form studied, to the second-order, was semi- 


circular. The accuracy of the second-order solutions was 


eel 








FIGURE 14. Wave Amplitude Generated by a 
Plunger Type Wave Generator 








FIGURE 15. Wave Amplitude Generated by a 
Flapper Type Wave Generator 
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established by comparing with the work of other researchers. 
Moemll7,10), Parissis [27], and Potash [31] have reported 
numerical solutions for second-order exciting force coeffi- 
cients , oc? and phase angle he), Figure 16 compares the 
values of ee and a. KepelEvcG Dye LeC.erarissis, and 
metasiewith the numerical results obtained by the FEM in this 
work. The points presented from the other researchers' works 
feeemeel Obtained from Potash [31]. Potash indicates Lee's 
points are corrected in some manner from his original work 
[17] and otherwise unpublished. 

Figure 16 shows excellent agreement for ay ancey cic, 
good agreement for ae Only Potash's points show signi- 
Mmecant deviation and then only for 56 near 1.5. Potash 
reported some numerical difficulties for a range of 6 near l 
and therefore Potash's results are not presented for those 
values of'6 which are suspected to be in error. Similarly., 
figure 17 shows a comparison of y 62) (phase angle) results. 
The figure shows the agreement is not as good. However, none 
of the other works shown (Lee, Parissis, Potash) agree in any 
consistent fashion over a large range of 6. The reason for 
the discrepancy has not been resolved. 

2. semi-Circular Hull — Finite Depth 

Figure 18 presents the effect of depth on pe and 

_? as obtained by the FEM in'’this work. The author knows 


of no other solutions. Solutions were obtained for values 


Sena trom 6.0 (infinite depth) to 1.5. The figure shows 


ie 













FIGURE 16. 


® Potash(3l) & 
we ere, ee 

-¢ Lee (17,|8) 

& Parissis (27) 


Theory 


second-Order Exciting Force Amplitude in 
heave moc orreular Hull, Infinite Depth 
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Theory 
e Potash(3l) 
-¢ Lee(I7,!8) 
yx Parissis(27) 





FIGURE 17. Second-Order Exciting Force Phase Angle 
in Heave, Semi-Circular Hull, Infinite Depth 
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FIGURE 18. Variation of Second-Order Exciting 
Horce Amplitudes with Depth, 
Semi-Circular Hull in Heave 
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that ? and a are most sensitive to finite depth at high 


values of 5 or in very shallow water (h/d < 1.5). 


Figure 19 shows the corresponding effect on y 62) out 


- 


finite depth for the same range of h/d. Note that at 
(2) 


~~ 


mor= 1.4 6 has negligible effect on y 


Also the effect of decreasing depth causes a monotonic 


(2) 24b - al seeveln®) - 


weelacion of yY 





FIGURE 19. Variation of Second-Order Exciting Force Phase 
Angles with Depth, Semi-Circular Hull Bhi. eleelnze 
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Migure secO presents the ratio of am to p(l) 
plotted vs 6 for various values of h/d. The curves peak 
close to 6 = 1 for all values of h/d studied. The curves 
peak because the first-order force Fi) passes through a 


minimum very near 6 = 1 for all values of h/d studied. 





FIGURE 20. Ratio of Second-Order Exciting Force Amplitude 
40 First-Order Exciting Force Amplitude, 
Semi-Circular Hull in Heave 
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A quantitative illustration of the physical signif- 
icance of the nonlinear effects (Second-order) is readily 
made using figure 20. Assume a ship with an approximately 
eeeeular hull form is being excited such that the maximum 
acceleration of the motion is 0.2 g's. The maximum time- 


dependent amplitude of the force is given by 


res eic| + en|e'*)| (130) 


which may be written 


a(2 
BC” nex =lte te | (131) 
e f F 
Assume h/d = 1.5 and 6 = 1.0, the maximum ship acceleration 


(in g's) may be written 


_ Se e 6. (132) 


E 


Equation 132 implies e« = 0.2 under the assumed conditions. 


Using figure 20 and equations 131 and 132 


Gt) | 


e f 1 


ot (Onl) Clea) = 1.35 (133) 


Therefore, neglecting the second-order contribution could 
result in a thirty-five percent error in |F(t)|_., 
Figure 21 further demonstrates the nonlinear (second- 


Seder) ettecets on F(t). The function R(t)/e £6}? woot bed 
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mec. The ship's displacement as a function of time is also 
shown on the plot for reference (cosine function of unit 
amplitude). The parameters chosen for the plot are h/d = 1.5, 
§6 = 1.0 and €« = 0.2 (corresponds to an equivalent acceleration 


mo .c g's). 


2.0 





€=0.2,0#1,0 
h/d 21.5 


10 2G  #-232 40 


MucCluREeol otal Exciting Force vs. Time 
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Figure 22 illustrates, for the same values of the 
Memeameters, the nonlinear effects on the gravity wave 
generated by the ship's disturbance. The figure presents 
the sum of the dimensionless first-order wave amplitude and 
the second-order Stokes' wave (both have celerity C1). 
Separately, the second-order wave from g(e) (celerity C5) is 


plotted. Both curves are plotted vs x/b for fixed time. 





FIGURE 22. Second-Order Wave Form Effects 


87 





VI. CONCLUSIONS 


It 1s concluded that this work has demonstrated the 
feasibility of using the finite element method to solve the 
mimear and nonlinear probiems of the type presented. Further, 
fmientinite element method is shown to be quite versatile in 
that various boundary conditions may be applied to the 
boundary value problems studied, with only minor changes in 
the computational scheme. In addition, the isoparametric 
ememtus chosen are capable of representing arbitrary hull 
Shapes. 

mes iurther concluded that A Pie OleGdew sso ululcy 7 Otenmlic 
effect of finite depth on the ship motions studied demon- 
Pemeaves that little change in the solutions occur. However, 
ae 1S asserted that the second-order results obtained 
demonstrate that nonlinear (second-order) effects on heaving 
motions are most significant in very shallow water; particu- 
darly for motions at a frequency corresponding very closely 
to the natural frequency of free oscillations of the floating 
body. 

ines present state of digital computer capacity appears to 
Preclude the extension of the present problem to three- 
dimensions. However, the extension of this study to other 
modes of motion (sway and roll) is possible as demonstrated 
Merecpne dinear sway problem. Further, the possibility of 
studying non-symmetric, two-dimensional hull forms is 


peecested. 
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An adequate finite element mesh to represent the region 
for the nonlinear problems studied tends to tax the core 
moeeeipy Of the computer used. Nevertheless, meaningful 


results were obtained through careful design of the computa- 


tional scheme. 
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APPENDIX A 


Complex Algebra 


The development which follows determines the correct 


form of the complex numbers Wy and Wo in the following 


expression 

Re{z i iRe(z eneiey = Re{w,. + w Bo +ots (Al) 
if °C al ° 

Z. and Zo are complex numbers with moduli eT and Py and 


arguments 6. and 05. It follows using appropriate trigono- 


1 


metric identies that 


On Vg. 
Retr e eo" }Re{roe aoe 


Me ‘ Oy a 
rie [cos6,cos6., let COS Ot met siné,sin6, Fle Oa) 


2 


rr 
; ; om ie 
- sin(6,+6,)4sin2ot] = 5 [cos(6,-8,) + cos (6,+65t+20t) J Chee) 





Equation A2 implies that 








Z4 Z 
—— , 0 Wy ee (A3) 


The symbol (~) implies conjugation. 


Equation A3 shows the origin of the factor % in various 


peaees in the text. 
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APPENDIX B 


Isoparametric Elements 





x,y plane 


FIGURE Bl. Element Mapping from €,n Plane to x,y Plane 


The essential idea in the isoparametric element lies 
in the following equations which define a map from the &,n 


Parene tO the x,y plane. 


x =N x, (B1) 


SNS (B2) 


The vector N© is a 1x12 row vector of element shape functions 


one for each node in figure Bl. The vectors x= and y< are 


Sal 





Hex. cOlumn vectors representing the respective x and y 
coordinates of the nodes in the x,y plane (figure Bl). 
The shape functions are conveniently defined using the 


Menventions of Zienkiewicz [45]. Let 


Ba as 5 lil, & ule , (B3) 


where € and n are the coordinates of any point of the square 
element in the €,n plane. bs and n, are the coordinates of 
the ue node. The shape functions are given by the following 


set of equations. 


Corner Nodes: 


NS = py (1 + E,)(1 # ng)E-10 + 9(E* + n®)] , (Bl) 
Edge Nodes: 
= _— up 

for €. = +1, ns te 

NG = = (1 + E)(2 - n°) (1 + 9g) (B5) 
1 : 

for §&, = #3, ny = Fl, 

NE = se (1 + ng) (2 ~ 6°) (1 + 964) (B6) 
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The element is defined to be isoparametric when the same 
shape functions are used to achieve the element mapping as 
Pemuisced CO define the field function in the x,y plane. 

There are certain consistency requirements Bn the element 
shape functions as discussed by Zienkiewicz [45]. The shape 
functions chosen satisfy them all. 

The mapping has the advantage of being able to represent 
a region with curvilinear sides. The area and line integrals 
given in part III are calculated by the expressions given 
below. 


From the calculus, the following relation is known to hold 


Ne . 

7 = J ~ OBC) 
e ae e 

N N 

~N) ~JY 


feercechne real 2x2 Jacobian matrix of the transformation. 


oer 


ee ee 
N N 
a oe piece 
J = (B8) 
: eae ee 
~~ ~N~ | 
it follows, provided J is not singular, that 
e e 
N 
N Pies 
=a (B9) 
ne _ ne 
~J ~T} 


a6 





Then, for example, equation 86 in part III may be written at 


the element level as 


e 
N 
e a el * 
H- = f Le Ny 1 dRe (B10) 
~ Re ~— . e 
N 
zai 


Meme equations B8, B9, and 
dRe = dxdy = |J| d&dn_ , (B11) 


equation B9 may be written 


LS: 
N 
aie il At 7 ate 
Ho = SY [NE Ne J Ca7*]* [a7] |J| agan . (B12) 
i -1 = = 
~1 


Equation Bl2 demonstrates the method of the evaluation 
of the matrix H© or in Scher die vine wmavit He Ihe integrations 


are accomplished uSing six-point Gauss quadrature. 





APPENDIX C 


CoOMmeucet eerocram Econom Zat ion 





PEEURE Cl. Schematic of Coefficient Matrix 


The economizations in the computer solution of the linear 


system of equations defined by equation 102 in part IV are 


achieved by observing certain properties of the matrix Koh). 


~e 


First all interior nodes of a given mesh contribute only real 


(1) 


numbers to K Liomtema lsocmeurie Clenodes On Sy. So and 3. 


(1) 


~s 


Only the nodes on Sy contribute complex entries in K 
Further only nodes on 53 and Sh contribute entries which are 
frequency dependent. 

The nodes on 2 3 and Sy myoLCcally comprise no more than 
Perc een percent of Kit), Further, the nodes on Sh alone 


~ 


ZO 





(1) 


comprise nor more than five percent of the total rows in K R 
Therefore if en? is rearranged with equations for nodes a 
Sy aeuethe bottom, equations for S 3 nodes directly above the 
Sy equations, and all other rows of the matrix filling the 
remaining part of ge) the resulting matrix would appear as 
schematically represented in figure Cl. The region designated 
with a Roman numeral I (uncrosshatched) indicates the portion 
of the matrix completely real and frequency independent. 
Region II indicates the portion of the matrix which has 
elements that are frequency dependent. And, Region III 
indicates the portion of the matrix which iS complex and 


(1) 


frequency dependent. If the matrix K is formed as shown 


in figure Cl, eighty-five percent (approximately) of it may 
be eliminated only once uSing real arithmetic. This is 
possible because the matrix ~5°Q, does not have to be added 
“to-Region II until after Region I has been reduced by Gauss 
elimination, and the same fact applies to the matrix 

(o/c, )D in Region III. Then, after Region I is eliminated, 
the matrix aah is added to Region II (after saving the 


(1) 


previous results) and K is reduced to upper triangular 


form through Region II using real arithmetic. Now, if this 
result is saved and then the matrix (io/c,)D is added to 
Region III only a small set of complex equations (usually 
mourn 25) is left to solve. 

Miperethertinal elimination, the resulting matrix is in 


Semoplere woper-triangular form and ready for a back- 


substitution process. The right-hand side vector is complex 
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and frequency dependent in general. However, the required 
imeormattOn to perform the forward substitution process on 


memo contained in the final form of x 6) (upper triangular). 
i (1) 


Miemenoare, the forward substitution on ibor 


(1) 


may be done 
after K is processed. Then a back-sSubstitution of the 
resulting system using mixed-mode computer arithmetic yields 
the final solution vector oo. 

A close examination of the process just described reveals 
that for a whole family of frequencies only a small portion 
of ay has to be reassembled and eliminated again. These 


Matrix properties contribute toward a sizeable saving in 


computer time. 


SH 





el. 


BIBLIOGRAPHY 


meamek, J. R., An Automatic Mesh Generator Using Two 


gee tnree-Dimensional Isoparametric Elements, M. 8S. 
Thesis, Naval Postgraduate School, Monterey, California, 


Ee 3 « 


mempousrd, Y., and Coubert, J.F., "Numerical Study of 
Maemotvory Or Non-Linear Waves," Proceedings of the 
International Symposium on Finite Element Methods in 


Hiow Problems, Swansea, January, 1974, pp. 283-288. 


Bai, K. J., "A Variational Method in Potential Flows 
with a Free Surface," University of California, Berkeley, 
College of Engineering, Report No. NA 7e-2, 1972. 


Swemauen, D. W., Ll, Motion of a Ship at the Free Surface, 
Mes. Thesis, Naval Postgraduate School, Monterey, 
Cmaiornia, 1970. 


Fontannet, P., "Theorie de la Generation de la Houle 
Cylindrique par un Batteur Plan (2© Order d'Approxima- 
fom) .," La Houillie Blanche, Vol. 16, 1961, pp. 179-197. 


imeem, W., "Oscillations of Cylinders in or Below the 
Free Surface," Naval Ship Research and Development 
Center, Report No. 2375, October, 1967, pp. 1-48. 


feummrioon, C. J., The Consistent Second-Order Theory 

of Wave/Structure Interaction, paper to be presented 
at the 14th International Conference on Coastal Engi- 
neering, Copenhagen, 24-28 June, 1974; To be published 
in the Proceedings of the above: September, 1974. 


femrrson, ©. J., "Hydrodynamics of Large Objects in the 
Sea Part I - Hydrodynamic Analysis," Journal of Hydro- 
Mmawrios, Vol. 8, No. 2, January, 1974, pp. 5-le. 


tame O., “Berechnung der durch Schwingungen Eines 
Schiffskorpers Erzeugten Hydrodynamischen Krafte," 
der Schiffbautech. Ges., Vol. 47, 1953, pp. 277-296. 


Havelock, T., Philosophical Magazine, Vol. 8, No. 51, 
October 1929, pp. 569-576. 


Havelock, T., "Waves Due to a Floating Sphere Making 


Peommodue Heaving Oscillations,” Proceedings Royal 
p@evwety of London, Vol. 231, 1955, pp. 1-7. 


98 





ee 


3 


14. 


Ie 


hoe 


iti. 


Lt 


i. 


208 


2 1 


Ca 


Ae 


folate. i., “Finite Hilements for the Computation of 
Hydrodynamic Mass," Proceedings of the Symposium on 


PiMiiverelement Lechnigues, Stuttgart , 1969, PeeeoaC9—5 41). 


motels Ol the Motion of Floating Bodies," 
Communications on Pure and Applied Mathematics, Vol. 3, 
ioe, pp. 45-101. 


Kim, C. H., "Hydrodynamic Forces and Moments for Heaving, 
Swaying, and Rolling Cylinders on Water of Finite Depth," 
goumnal of Ship Research, Vol. 13, No. 2, 1969, 

pee l5/-154. 


Lamb, H., Hydrodynamics, Cambridge University Press, 
Sem edition, 1932. 


Landweber, L., and Macagno, M., "Added Mass of a Three 
Parameter Family of Two Dimensional Forms," Journal 
of Ship Research, Vol. 2 , No. 4, 1959, pp. 36-48, 


Lee, C. M., "The Second-Order Theory of Cylinders 
Secillating in a Free Surface," University of California, 
Berkeley, College of Engineering, Report No. NA-66-7, 
September, 1966. 


Lee, C. M., "The Second-Order Theory of Heaving Cylinders 
in a Free Surface," Journal of Ship Research, Vol. 12, 
mIOC, Pp. 313-327. 


Lewis, F. M., "The Inertia of Water Surrounding a 
feeerating Ship,” Transactions of the Soctety vofteNavall 
Architects and Marine Engineers, Vol. 37, 1929, 

pp. 1-20. 


MeeCamy, R. C., “A Source Solution for Short-Crested 
Waves," La Houille Blanche, No. 3, July-August, 1957, 
Peet, 319-509. 


Matsumoto, K., “Application of Finite Element Method 
mommogded Virtual Mass of Ship Hull Vibration,” Journal 
mene DOCLeLTY Of Naval Architects of Japan, Vol. le/7. 
eee. 03-90. 


Matsuura, Y., and Kawakami, M., "Calculation of Added 
ercdeal Mass Moment of Inertia of Ship Hull Vibration 
by the Finite Element Method," Journal of the Society 
GmeNiaval Architects of Japan, Vol. 124, MOGI. Op: 2ol- 


291. 


Neovmeouem ire to, Chedauly, Dp. W., LL, and Smith, D. A., 
Jr., "Finite Element Solution for Added Mass and Damping," 
Proceedings of the International Symposium on Finite 
Pelemene Methods in Flow Problems, Swansea, January, 


1974, pp. 159-170. 


Si, 





a4, 


256 


26. 


21 
ZO, 
2 


SUF 


Be. 


Bc. 


Se 
34. 


BD. 


Newton, R. &., Finite Element Analysis of Two-Dimensional 
maged Mass and Damping, to be published in Collection 

of Papers presented at the International Symposium 

on Finite Element Methods in Flow Problems, Swansea, 

J. Wiley and Sons, London. 


Ogilvie, T. F., "First and Second-Order Forces on a 
Cylinder Submerged Under a Free Surface," Journal of 
Piad Mechanics, Vol. 16, 1963, pp. 451-472. 


Paee we lets, and Tuck, EB. 0., "A Rational Strip 
micery ef Ship Motions,” University of Michigan, 
Department of Naval Architecture, Report No. 013, 1969. 


Parissis, G., "Second-Order Potentials and Forces for 
Oscillating Cylinders on a Free Surface," Massachusetts 
Institute of Technology, Report No. 66-10, 1966. 


Pauline, J. R., and Richardson, R., “Measurement of 
Pressures, Forces, and Radiating Waves for Cylinders 
Oscillating in a Free Surface," University of California, 
Berkeley, Report No. 82-28, 1962. 


mauling, J. R., and Porter, W. R., “Analysis and 
Measurement of Pressure and Force on Heaving Cylinders 


jm a Free Surface," Proceedings of the Fourth National 
Womeress Of Applied Mechanics, 1962, pp. 1309-l3060e. 
Porter, W. R., "Pressure Distributions, Added-Mass, 

and Damping Coefficients torr Gy Linders Oscillating in a 


' Free Surface," University of California, Berkeley, 


Report No. 83. Nene G0] 


Potash, R. L., "Second-Order Theory of Oscillating 
Peaiders, University of California, Berkeley, Report 
No. NA70-3, 1970. 


Remenee. M. @., Llmpact of Pinite Element Techniques 
Pome ractical Design of Ships Structures," Proceedings 
of the Symposium on Finite Element Techniques, Stuttgardt, 
Oe eo wm 09-551. 


Salveson, N., On Second-Order Wave Theory for Submerged 
ie Dimenatona. Bodies, Phe D. Dissertation, University 
of Michigan, Ann Arbor, 1966. 


Poumeers, tl. &., Hydrodymamics in Ship Design, Vol. ed, 
The Society of Naval Architects and Marine Engineers, 
ieumeuwork, 195/. 


weer. oe Water Waves, Interscience Publishers, 
Inc., New York, 1957, pp. 19-54. 


100 





36. 


37. 


38. 


Boi. 


HO, 


Al. 


42, 


USe 


Hy, 


Me 


H6, 


mood, fs, “Damping Forces and Added Mass of Ships 
Heaving and Pitching," Kyushu University, Research 
Institute of Applied Mechanics, Report, Vol. 8 
eesti, 1959; Vol. o, pp. 36-39, 1960. 


3 


Tuck, E. 0O., "The Effect of Non-Linearity at the Free 
Surface on Flow Past a Submerged Cylinder," Journal 
of Fluid Mechanics, Vol. 22, 1965, pp. 401-414, 


Meee ll, K., “On the Heaving Motion of a Circular Cylinder 
on the Surface of a Fluid," Quarterly Journal of Applied 
Mechanics and Applied Mathematics, Vol. 2, No. 2, 1949, 
Preece o-231. a 


Messer, W., and van der Wilt, M., "A Numerical Approach 
Memene ovudy of Erregular Ship Motions," Proceedings 

Or the International Symposium on Finite Element 
Methods in Flow Problems, Swansea, January, 1974, 

hee 277-202. 


Vugts, J., "The Hydrodynamic Coefficients for Swaying, 
meee wie, and Rolling Cylinders in a Free Surface," 
Memecrlands ohip Research Center, Report No. 1125, 
1968. 


Wehausen , J.V., and Laitone, Ee VerecouictacesNaves 
mieyelopedia of Physics, Vol. 9, edited by S. Flugge, 
Springer-Verlag, Berlin, 1960. 


(ieee. , and Ursell, F., “Surface Waves Generated 
peeeom Oscillating Cylinder on Water of Finite Depth: 
Theory and Experiment," Journal of Fluid Mechanics, 
femeell, No. 4, 1961, pp. 529-551. ; 


Zienkiewicz, O. C., and Newton, R. E., "Coupled Vibra- 
tions of a Structure Submerged in a Compressible Fluid," 
Proceedings of the Symposium on Finite Element Techniques, 
Biameueort, 1969, pp: 359-379. 


Ptemiciewilcz, O. €., Irons, B. M.; and Nath, B., "Natural 
Frequencies of Complex Free or Submerged Structures by 
the Finite Element Method," Proceedings of the Conference 
Semeesoreacwons in Civil Engineering, Institution of 


cweeneineers, London, 1965. 


weewierewicz, O. C., The Finite Element Method in 
Fngineering Science, McGraw-Hill, London, 19/1. 


miemicrewie,, O. C., and Cheung, Y. K., "Finite Elements 


Meeonemoolution of Field Problems,” The Engineer, 
September, 1965, pp. 507-510. 


LOL 





eNom oukLON Liow 


No. 


Defense Documentation Center 
Cameron Station 
Pexandria, Virginia 22314 


Paprary, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Professor R. E. Newton, Code 59Ne 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


Commander, Naval Shipyard. 
Mie: LT David A. Smith, Jr., Code 338.22 
bong Beach, California 91/701 


miaoressor J. R. Paulling 
Department of Naval Architecture 
University of California 
Berkeley, California 94720 


Naval Ship Research and Development Center 


Engineering Research Project in Hydromechanics 


Bethesda, Maryland 20034 


Naval Ship Systems Command Headquarters 
Code 03 

Navy Department 

Washington, D. C. 20360 


OZ 


Copies 


it 
























| a ene 
anecic sults 


$9759 Smith 
€ wil Finite element anal- 


ysis of the forced os-~- 
cillations of ship 
hull forms. 





thesS5759 
Finite element analysis of the forced os 


8 
L DUDLEY KNOX LIBRARY 


